Mott metal-insulator transition on compressible lattices 
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The critical properties of the finite temperature Mott endpoint are drastically altered by a coupling 
to crystal elasticity, i.e., whenever it is amenable to pressure tuning. Similar as for critical piezo- 
electric ferroelectrics, the Ising criticality of the electronic system is preempted by an isostructural 
instability, and long-range shear forces suppress microscopic fluctuations. As a result, the endpoint 
is governed by Landau criticality. Its hallmark is thus a breakdown of Hooke's law of elasticity with 
a non-linear strain-stress relation characterized by a mean-field exponent. Based on a quantitative 
estimate, we predict critical elasticity to dominate the temperature range AT*/T C — 8% close to 
the Mott endpoint of k-(BEDT-TTF) 2 X. 

PACS numbers: 



Strong repulsion between electrons in a solid contain- 
ing approximately one electron per lattice site promotes 
insulating behavior as the electrons' motion is inhibited 
by the large energetic cost of having a site doubly occu- 
pied. A so-called Mott insulator is favored if the on-site 
Coulomb repulsion U exceeds the kinetic energy W while 
metallic behavior prevails for U/W <C 1, allowing for a 
first-order metal-insulator transition at a critical ratio of 
U/W. Usually, this critical ratio can be controlled by 
an external control parameter like pressure, p, or doping, 
giving rise to a line of first-order transitions in the (p, T) 
phase diagram plane, where T is temperature. This line 
of transitions terminates in a second-order critical end- 
point at a finite temperature T c beyond which the system 
can be smoothly transformed from the insulating to the 
metallic regime by varying p and T, see Fig. [T] 

The nature of the finite- T critical endpoint of the Mott 
transition attracted some attention recently. From gen- 
eral considerations, one expects it to belong to the Ising 
universality class 0, 0] similarly to the endpoint of the 
liquid-gas transition. The double occupancy of a single 
site plays here the role of the local Ising order param- 
eter of the transition. Measurements of the electrical 
conductivity on Cr-doped V2O3 [3| have confirmed this 
expectation. However, transport [H and NMR || mea- 
surements on the quasi two-dimensional organic charge- 
transfer salt k-(BEDT-TTF) 2 X JrJ-Q questioned this in- 
terpretation and suggested a different universality class. 
Subsequently, various theories were proposed to account 
for the unconventional behavior [1CH13I|. In particular, it 
was demonstrated [l4[ that the analysis of the conductiv- 
ity is intricate as its scaling exponents are not necessarily 
directly related to the scaling dimension of the Ising order 
parameter. Taking this into account, the conductivity 
experiments could be reconciled with Ising universality. 

In contrast to transport quantities, thermodynamics 
in principle allows for a straightforward interpretation 
in terms of a standard critical scaling analysis. Ultra- 
sound experiments on k-(BEDT-TTF)2X revealed a pro- 
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FIG. 1: The line of first-order Mott transitions termi- 
nates at a finite temperature critical endpoint. The pressure- 
tuned endpoint exhibits Landau criticality due to the non- 
perturbative Mott-elastic coupling (yellow regime). The 
dashed lines bound possible low temperature phases like an- 
tiferromagnetism and superconductivity. 



nounced softening close to the Mott endpoint [8j , but an 
analysis of the critical behavior has not been performed 
yet. Thermal expansion measurements 15|, [16| on the 
same system are consistent with two-dimensional critical 
Ising behavior even though an experimental verification 
of scaling exponents was not possible until now. 

While the Mott transition of an idealized correlated 
electron system should exhibit Ising criticality, a cou- 
pling of electrons to crystal elasticity drastically changes 
its critical properties. The Mott transition is very sensi- 
tive to the presence of elastic strain in the atomic crystal 
lattice as it alters the overlap integrals of electron wave- 
functions between adjacent lattice sites. This in turn 
changes the kinetic energy W, allowing for an efficient 
tuning of the transition by applying external stress, e.g., 
a compressive pressure as in Fig. [TJ Conversely, the criti- 
cal Mott system exerts an internal pressure on the elastic 
system to which the crystal lattice responds. The detec- 
tion of this response with the help of dilatometric mea- 
surements, e.g., thermal expansion, is a convenient and 
sensitive probe of Mott criticality. Sufficiently far away 
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from criticality, this lattice response is perturbative and 
the critical behavior itself remains unaffected. 

However, close to the Mott endpoint, the lattice nec- 
essarily reacts in a non-perturbative manner to the in- 
ternal stress, leading to a vanishing elastic modulus and 
thus to a breakdown of Hooke's law of elasticity. This 
was noted before by Krishnamurthy and collaborators 



171 1 1 81 ) within the framework of the compressible Hub- 
bard model. Importantly, we point out here that this 
breakdown of Hooke's law is generically accompanied by 
a crossover from Ising criticality to Landau critical be- 
havior with mean-field exponents. At the origin of this 
change of universality class are the long-ranged shear 
forces of the atomic crystal lattice. They become in- 
strumental as an elastic modulus becomes small so that 
eventually Landau mean-field behavior prevails close to 
the Mott endpoint. 

Within an effective field theoretic description, we con- 
sider a coupling of the elastic strain tensor Ey to the Ising 
order parameter <j> of the Mott transition, 



Ant = -ll,ij£i]4>+ ^2,ij£%j(f> 2 



(1) 



where 71 ij and 72, ij are elastic coupling tensors. Inter- 
estingly, as the Ising symmetry of the Mott endpoint is 
an emergent symmetry, a linear coupling Ji,ij of strain to 
the order parameter is generally allowed. The quadratic 
coupling j2,ij is less important, but for completeness we 
include it in the following discussion. The effect of a lin- 
ear coupling of an order parameter to strain was consid- 
ered by Levanyuk and Sobyanin [l9j | and independently 
by Villain 20] in the context of critical ferroelectrics, who 
showed that it suppresses critical long-wavelength fluc- 
tuations. As a consequence, for sufficiently large 7i,y, 
the Ginzburg criterion is never fullfilled, thus stabilizing 
Landau mean-field behavior. 

An analysis of the effective elastic Hamiltonian [21. 22| 
suggests that even if the critical subsystem is controlled 
by an interacting renormalization group fixed point, a 
small linear elastic coupling 71,^ can recover mean- field 
behavior sufficiently close to the transition. It turns out 
that the singularities associated with the Mott endpoint 
induce via the coupling ji^j a macroscopic instability 
of the crystal lattice. At such a lattice instability, an 
elastic modulus associated with the macroscopic strain 



solid. The phonon velocities are determined by the 3x3 
matrix Mjfe(q) = J^ji Cijkiqjqi depending on momentum 
q, and its eigenvalues generally remain positive even if 
an eigenvalue of Cijki vanishes. Possible exceptions may 
be acoustic modes with momenta in certain lattice di- 
rections. For general momenta, however, the acoustic 
modes remain non-critical and, as a consequence, the 
structural transition is described by Landau's mean-field 
theory [5l|,[22|. 

Neglecting these non-critical acoustic modes, the 
macroscopic strain Eij is determined by the effective po- 
tential 
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where Cyjy is the elastic constant matrix in the absence 
of Mott-elastic couplings 7 n ,jj 5 with n — 1,2, and (jjj is 
an externally applied macroscopic stress. The free energy 
density /sing is attributed to the critical electronic subsys- 
tem and is governed by the Ising universality class. The 
two relevant perturbations ho and to quantify the dis- 
tance to criticality (for J n ,ij = 0) and generally depend 
on temperature T . In order not to distract with cumber- 
some notation and to focus on the mechanism at play, let 
us assume that the electronic subsystem mainly couples 
to a certain singlet, E, of the irreducible representations 
of the crystal group. We can then limit ourselves to an 
effective potential for E only, 



V(E) = ^ 



£p+/sing(t0+72#,/l0+7l£)- (4) 



Here Ko is the corresponding modulus for j n = 0. Fur- 
thermore, we considered for simplicity the application of 
a hydrostatic pressure, <Jij = —pSij, assuming a finite 
overlap with the singlet E. The thermodynamic free en- 
ergy density obtains after minimizing this potential with 
respect to E. 

The sensitivity of the Mott endpoint with respect to 
pressure tuning becomes manifest if the potential ((3]) is 
minimized perturbatively in 7„. In zeroth order one has 
E = p/Ko so that the free energy density becomes 

P 2 

Fpert = + /smg(*0 + l2p/K , h Q + Jip/K ). (5) 



matrix, C D 



the elastic strain also contains a part that carries finite 
momentum and describes the long-wavelength acoustic 
modes 



vanishes 1231. In addition to Ei 
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The elastic coupling induces a pressure dependence of the 
field E^, i.e., an eigenvalue of the 6 x 6 elastic constant arguments of the function / sing which allows to control 

the distance to criticality by varying p, thus enabling 
pressure-tuning of the Mott transition. 

However, it is important to realize that such a per- 
turbative treatment necessarily breaks down sufficiently 
close to the endpoint. This becomes evident after ex- 
panding the potential in a Taylor series, 

V(E) = f ains (t,h) 

-(p-p)5E+^-SE 2 + ^5E A + 0(5E 5 ), (6) 



£ij{v) = E^ + eij(r), 



(2) 



with J d 3 reij(r) = 0. Importantly, at a lattice insta- 
bility, the velocities of the acoustic modes soften, but 
generally remain finite due to the shear stiffness of the 
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FIG. 2: Phase diagram close to the finite- T Mott endpoint in 
the (a) absence and (6) presence of a linear elastic coupling 
71 7^ (and 72 = 0) ; to and ho are the two relevant fields and 
p is the pressure. Whereas for (a) the endpoint is Ising critical 
for (b) a crossover is induced from Ising to Landau mean-field 
criticality for \to — t c o\ < t c o and \p — p c \ < Ap*. In the 
uncolored and yellow shaded regime in (6) thermodynamics 
is governed by Eq. ((5JI and Eq. JB), respectively, and in the 
brown shaded regime the full potential Eq. @ must be used. 



where 5E = E-E,t = t + j 2 E and h = h + 71 The 
value of E is conveniently chosen such that the prefactor 
of the cubic term, SE 3 , in the expansion just vanishes. 
The pressure p reads p = K E + (726% + 7i9^)/ s i ng (f, h), 
and the quartic coupling u is given by fourth-order 
derivatives of / S i ng - Importantly, the modulus gets renor- 
malized by the susceptibilities \ab = —9 a dbf S ing{i, h) 
with a, b = t, h, 



K = K Q - ^Xhh - 27i72Xfcf 



' llxtt- 



(7) 



The most singular susceptibility is Xhh which necessarily 
diverges, XhK 00 ; as the endpoint (t, h) = is ap- 
proached. Hence, irrespective of the magnitude of the 
linear elastic coupling, 71 7^ 0, the divergence of Xhh wm 
drive the effective modulus to zero at a finite value of 
i where the Taylor expansion in Eq. ([SJ is well defined. 
The resulting isostructural instability at K = and p = p 
identifies a mean-field endpoint in the phase-diagram. 
This critical endpoint at (t c o,p c ) preempts the Mott-Ising 
endpoint in Fig. [2jb). The coupled Mott-elastic system 
thus avoids the Ising singularities by developing a non- 
perturbative strain response to pressure changes. Mini- 
mization of Eq. (JB]) for K — yields SE — (6(p—p)/u) 1 / s 
with the Landau value 6 = 3, clearly violating Hooke's 

law of elasticity. This violation sets in for \p — p\ < Ap* 

3/2 



at K = with Ap* = Kg' V 6 / M - 

As a concrete example, we assume that the critical 
electronic subsystem is effectively two-dimensional, and 
the function / s i ng in Eq. (j4j is determined by the 2d Ising 
model 

Wa) = /o ^log^ + ifti 18 / 16 *^^!- 8 / 16 )) , (8) 

where /o has the dimension of a free energy density (and 
t and h are assumed to be dimensionless) . Using the re- 
sults of Ref. [24| , the scaling function $ can be evaluated 




FIG. 3: The second derivative —d1 T as a function of to ex- 
hibits a pure mean-field jump at criticality due to the linear 
Mott-elastic coupling (solid curve) that preempts the Ising 
singularity (dashed curve), see text. 



numerically. In order to illustrate the mean-field charac- 
ter of the shifted Mott endpoint, we show in Fig. [3] the 
second derivative —d^T of the free energy density T for 
71 7^ (and 72 =0). The pressure is fixed to the critical 
value p c so that the endpoint is crossed as a function of 
to, i.e., along the vertical axis in Fig. [U[b). The solid 
curve shows the behavior obtained from minimizing the 
full potential ((4]) while the dashed curve follows from the 
perturbative expression J5]) . The latter exhibits the char- 
acteristic logarithmic divergence of the 2d Ising model at 
to = 0. However, the non-perturbative renormalization 
of the elastic constant results in a preemptive mean-field 
transition at to c > so that the logarithmic divergence is 
cut off and — df iF instead shows a mean-field jump and 
remains finite [19j , 

These considerations are directly relevant for k- 
(BEDT-TTF) 2 X close to its Mott endpoint. In Ref. [ijj, 
the perturbative free energy density ([5]) was used to- 
gether with Eq. for the interpretation of thermal ex- 
pansion measurements. For the so-called <i8-Br crystal 
#1 in Ref. [3], for which p — p c ~ 50 bar at ambient 
pressure and T c ks 30 K, the following fitting parameters 
were obtained: f ~ 5.7 bar, ho + ^ip c /K ~ -0.004(T- 
T c )/T c , and 71/-K0 ~ 0.07/ kbar, where the scaling free- 
dom was exploited to choose to + l2Pc/ Kq = (T — T c )/T c . 
As the exact critical temperature T c of the crystal is not 
known, there is no reliable estimate for 72. In the fol- 
lowing, we neglect the subleading corrections due to 72 
and use 72 = 0. A crucial question concerns the ex- 
tension of the non-perturbative regime in order to assess 
whether an experimental investigation of the crossover 
from Ising to Landau criticality is feasible. With the 
above fitting parameters and the estimate for the bare 
modulus Ko ~ 122 kbar we can estimate the width 
of the Landau critical regime in pressure, Ap*, and tem- 
perature, AT* = t c oT c for this compound, see Fig. Q] 



Ap*«45bar AT*«2.5K. 



(9) 
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FIG. 4: Lattice strain SE as a function of applied pressure 
close to the critical temperature T c . At T c (red line), the 
strain is linear in the applied pressure for \p — p c \ > Ap* 
(dashed-dotted line asymptote) but becomes non-linear close 
to the endpoint, where SE ~ \p — p c | 1/ '' 5 with the mean-field 
exponent 8 = 3 (dashed line asymptote), signaling a break- 
down of Hooke's law. 



These values are sufficiently large to allow for an exper- 
imental detection of the crossover phenomena. In fact, 
the d8-Br crystal seems to be located already within the 
crossover regime as Ap* is on the same order as the dis- 
tance p — p c [Hf . 

With the above fitting values and the value for Kq 
we can predict the thermodynamics with the help of the 
potential (gj) and Eq. ©. In particular, the crossover 
at Ap* is illustrated in Fig. |4] which shows the expected 
lattice strain as a function of applied pressure. Far away 
from the transition, \p— p c \ 3> Ap*, the strain is linear in 
the applied pressure, thus obeying Hooke's law. However, 
at T c the pressure-strain relation becomes non-linear for 
p — > p c with mean-field exponent 8 = 3. This breakdown 
of Hooke's law and the concomitant divergence of the 
associated modulus serves as a smoking-gun criterion for 
the detection of the Landau critical regime where the 
Mott-elastic coupling becomes non-perturbative. 

In Fig. EJa) we show the thermal expansion, a — 
dpdrJ 7 , as a function of temperature for different pres- 
sure values (solid lines) . For comparison, the dashed lines 
demonstrate the corresponding Ising critical behavior ob- 
tained from the perturbative expression of Eq. ([5]). The 
latter is a good approximation far away from the end- 
point but fails close to it and, in particular, exhibits a 
peak at a temperature ~ T c — AT* , that is smaller than 
T c , see also Fig. [2] The crossover is identified when the 
solid and dashed curves at a given pressure start to devi- 
ate substantially. Finally, Fig. [5jb) displays the thermal 
expansion as a function of pressure for different temper- 
atures. Note that the sign change of the thermal expan- 
sion in Fig. [5jb) can be related to entropy accumulation 
similarly as in the case of quantum criticality [jfjj ]. 

In Ref. [3] the conductivity, a, was interpreted to scale 
with the energy-density of the Ising model, a ~ dt fa- 
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FIG. 5: Singular part of thermal expansion (solid line) as a 
function of (a) temperature for different pressure values and 
of (b) pressure for different temperatures using the estimate 
Tc w 30 K. The perturbative behavior (dashed line) is a good 
approximation away from the endpoint, \p — p c \ 3> Ap* and 
IT - Tel > AT*. 



If this interpretation holds across the crossover discussed 
here, i.e., a ~ d^T, one would also expect signatures in 
transport at scales AT* and Ap* of Eq. ©. Interest- 
ingly, whereas the pressure dependence of a measured 
in Ref. [4j does not show such a signature, there are 
indications of a crossover in cr(T) at around T c ± 1 K. 
Clearly, detailed dilatometric studies are favored to iden- 
tify unambiguously the crossover to Mott-Landau criti- 
cality. Our estimate for k-(BEDT-TTF) 2 X, Eq. ©, indi- 
cates that this is experimentally feasible. This identifies 
this compound as a promising candidate to investigate 
the strong coupling between electronic and elastic de- 
grees of freedom close to the finite-T Mott endpoint and 
the concomitant change in universality class. 
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Paul and A. Rosch. This work is supported by the DFG 
grants SFB 608 and FOR 960. 
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